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The case of ac transport at in-phase alternating applied magnetic fields for a superconduct- 
ing rectangular strip with finite thickness has been investigated. The applied magnetic field 
is considered perpendicular to the current flow. We present numerical calculations assuming 
the critical state model of the current distribution and ac loss for various values of aspect 
ratio, transport current and applied field amplitude. A rich phenomenology is obtained due 
to the metastable nature of the critical state. We perform a detailed comparison with the 
analytical limits and we discuss their applicability for the actual geometry of superconduct- 
ing conductors. We also define a loss factor which allow a more detailed analysis of the ac 
behavior than the ac loss. Finally, we compare the calculations with experiments, showing 
a significant qualitative and quantitative agreement without any fitting parameter. 

I. INTRODUCTION 



The behavior of a superconductor transporting an alternating current or exposed to a magnetic 
field varying in time have been a wide subject of study from the early 1960's jl|, . However, 

the case of a simultaneous alternating transport current and applied magnetic field remains unclear. 
This situation is found in superconductor windings where each turn feels the magnetic field of all 
the other. Windings are present in many applications, like ac magnets, transformers and motors 
0,0,12,0]. From the practical point of view, it is of fundamental importance to understand, predict 
and, eventually, reduce the energy loss (or ac loss) in the superconductor. Actually, the reduction 
of the ac loss is vital for the applicability of superconductor electrical technology [3, O]. The 
study of the ac loss is also interesting for material science, since in can be used for characterizing 



superconducting samples 
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The superconductors suitable for electrical applications are hard type II ones Jdl- Nowadays 
there is a great scientific effort in the development of silver sheathed Bi2Sr2Ca2Cu30io (Ag/Bi- 
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2223) tapes and YBa2Cu307^ (YBCO) coated conductors, which are high-temperature super- 
conductors, and MgB2 wires These superconducting tapes and wires have a cross-section 
roughly rectangular or elliptical. In this work we will consider wires with rectangular cross-section 
(or rectangular bars), leaving those with elliptical cross-section for further studies. We also restrict 
to the situation when the ac applied field is uniform and in phase with the transport current. 

Hard type II superconductors can be well described by the critical- state model (CSM) proposed 
by Bean and London Pi , which assumes that the magnitude of the local current density cannot 
be higher than a certain critical value Jc. 

For the situation of only transport current, the CSM was first applied by London and Hancox 
in the early 1960's in order to analytically describe simple geometries, like infinite cylinders and 
slabs, P, Later, an important step forward was done by Norris, who analytically deduced the 
current distribution and the ac loss for an infinitely thin strip by means of conformal mapping 
transformations ll™ . The case of a strip with finite thickness can only be solved numerically, as 
done by Norris [isj l. Fukunaga et al ^3]) Daumling Q| and Pardo et al 1^ . 

The first analysis about the CSM with only ac applied field was done by Bean for an slab with 



Ic |. The case of a thin strip with perpendicular applied field was 
20| following the Norris' technique ll5|l. The current distribution 



applied field parallel to the surface 

analytically solved by Brandt et al ^ ^ ^ ^ 

for a strip with finite thickness were numerically calculated by Brandt [21( and Prigozhin [2^ > and 
the ac loss by Pardo et al Q]. 

Concerning the case of simultaneous alternating transport current and applied field, the most 
significant published calculations within the CSM are the following. In late 1970's Carr analytically 
derived the ac loss for an infinite slab in parallel applied field In the 1990's Brandt and Zeldov 
et al analytically calculated the current distribution in a thin strip using conformal transformations 
for the situation that the transport current and the applied field increase monotonically |2a, l2al. 
Al.ho... ..e. wo... p.ov.e ...en. f„...,ae. ..e. »e ac.* e„ M- mU, 
Brandt studies the values of transport current and applied field for which these formulae are valid, 
obtaining that they are not applicable for high fields and low currents ,2^. From that current 
distribution, Schonborg analytically calculated the ac loss for a thin strip |23]. 

For strips with finite thickness, there are no published works dealing with the simultaneous 
application of alternating currents and magnetic fields in a superconductor in the CSM, according 
to our knowledge. However, there are several theoretical works assuming a relation between the 
electrical field E and the current density J as E(J) = E'cd J|/ Jc)'^J/| J|, where Ec is an arbitrary 
value and n is a positive exponent. The current distribution and the ac loss are calculated in Refs. 
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and Refs. 



28 



2c 



31 



33l |. respectively, for several values of the alternating 



transport current and applied field. These published results are incomplete, not covering the whole 



range of combinations of ac current and ac fie^ 



that has been done for Ag/Bi-2223 tapes 



d. This contrasts with the extensive experimental 
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38i | and YBCO coated conductors 
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The objective of this paper is to provide a systematic study of the current distribution and ac loss 
for a superconducting strip of finite thickness assuming the CSM under simultaneous application 
of an alternating transport current and field. The effect of mainly three factors are considered: the 
aspect ratio of the cross-section and the aplitudes of the transport current the applied magnetic 
field. We also study the applicability of the CSM to actual superconducting tapes and wires of 
such shape. 

This paper is structured as follows. In Sec. ^ we present the numerical method used for 
the calculations and we discuss some general features, mainly about the energy of the system. 
The results and their discussion are presented in Sec. Ill II In Sec. IIVI the comparison with data 
measured in a high-temperature superconducting tape is reported. Finally, in Sec. 0we present 
our conclusions. 



II. NUMERICAL METHOD AND GENERAL CONSIDERATIONS 

Let us consider an infinitely long superconductor along the z axis with rectangular cross-section 
with dimensions 2a x 2b in the x and y directions, respectively, Fig. ^ The origin of coordinates 
is taken in the center of the strip. We study here the situation that the superconductor carries a 
sinusoidal time- varying current I{t) = cos simultaneously immersed in a uniform in-phase ac 
applied field Ha(t) = Hmcoswt oriented in the y direction. It is shown below fSecs. ITlDl and Hi El 
) that our results are not only independent on w but also on the specific time waveform of / and 
Ha, similarly to the case of only transport current or magnetic field 

In our calculations we will consider that first / and Ha are increased from zero to their maximum, 
starting from the zero-field cooled state of the superconductor. We call this process the initial stage. 
Following this stage, we regard the reversal one for which the current and applied field are decreased 
from Hm and 1^, respectively, to — Hm and — /m- Next, the applied field and current are increased 
back to their maximum, closing the ac cycle. We refer to this latter stage as the returning one. 
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A. The critical-state model in strips 

We assume that the superconductor obeys the CSM with constant critical-current density 
P]. The CSM corresponds to supposing a multivalued relation of electrical field E against current 
density J such that E = £'(|J|)J/|J| with an -BdJ]) that only takes finite values for \ J\ = Jc, being 
zero for \ J\ < Jc and infinity for |J| > Jc For an infinitely long strip along the z direction, 

the current density and the electrical field inside the superconductor are also in the z direction and 
they can be considered as the scalar quantities J and E, respectively. Although in principle Jl in 
the CSM can be lower than Jc, in a superconducting strip J only takes the values or itJc 42| . 

Let us start with the introduction of the main features of the current distribution in the initial 
stage for the case of transport current only (i.e. Jfa = 0) or when solely the magnetic field is 
applied (1 = 0). 

The behavior of a superconducting strip in the critical state model with Hg^ = and uniform 
Jc is detailed in In the initial stage, for any I > lower than the critical current, 

Ic = 4a6Jc, there exist a zone with J = surrounded by another one with J = Jc. The region with 
J = is usually called the current-free core. In this zone the electrical field is zero because in the 
CSM E{J = 0) = 0. With increasing I, the region with J = Jc monotonically penetrates from the 
whole surface inwards and the current-free core shrinks, until it disappears when / reaches Ic- In 
the CSM, / cannot overcome Jc since it is assumed that | J| < Jc. 

The situation when a magnetic fie 



any net current is described in 



21 



d is applied to a superconducting strip that is not transporting 
' |. For this case, the current distribution is antisymmetric 
to the yz plane. In the initial stage with i/a > 0, there are a zone with J = Jc in the right half 
and another one with J = — Jc in the left half expanding from the surface to a current-free core 
between them. Throughout this paper, we call the border between regions with different J as a 
current front. It is important to notice that in the current fronts J vanishes and, then, so does E. 
With increasing Hq^, the cross-section of the current-free core shrinks until it becomes a point at 
{x,y) = (0,0) at the characteristic field Hp, that is called the penetration field. At fields higher 
than Hp, the current distribution is the same as for Hg^ = Hp. 

When we now consider the simultaneous action of an applied magnetic field on a conductor 
transporting a nonzero current, one can expect that the qualitative behavior of the current distri- 
bution is similar to that one for only transport current or applied field. However, for some situations 
of I and J^a the current distribution presents a different behavior. We discuss this aspect in more 
detail below (Sec. IIIIA|) . 
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B. Minimization principle for the critical state model 



As discussed by several authors, such as Prigozhin 



et al ,45], and Sanchez et al |4(tI . 
assuming the critica 
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4ll |. Badia and Lopez 
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4^ 



, Bhagwat 



introduced in 
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41 



48| |. the distribution of current density for a superconductor 
state model is such that it minimizes a certain functional. The functionals 

lile in j^] it is proposed the magnetic energy as the 



44 



45l | are equivalent, w. 



quantity to be minimized. As shown in 
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4ll |. the principle of minimization of the functional, J-", 
can be derived from fundamental considerations. In Sec. Ill Cl we demonstrate that the minimization 
of is equivalent to minimizing the magnetic energy provided that in the initial stage the current 
front penetrates monotonically from the surface inwards. Some of the situations presented in this 
paper do not satisfy this condition; therefore we use the minimization of T as follows. 

Let consider the case of an infinitely long superconductor extended along the z direction carrying 
a transport current / and immersed in a uniform applied field in the y direction Hg^, Fig. ^ With 
this geometry, the current density is in the z direction and, therefore, so is the vector potential A 
if we assume the gauge V • A = 0. Then, we can regard these quantities as scalar. Following the 
notation of Prigozhin j4l| . the current at a certain time distributes in such a way that it minimizes 
the functional 



y^J{r)Aj{r)dS 



J(r)Aj(r)d5 



with the constrains 



+ / J(r)[Aa(r)-^a(r)]d5, 
Js 



J{r)dS 



\J\ < Jc 



(1) 

(2) 
(3) 



where S is the superconductor cross-section, Aj is the vector potential created by J, A^^ is the 
vector potential from the external field, and the quantities with hat correspond to those at the 
previous time layer. For infinitely long geometry, Aj can be calculated from 



^ Js'^ (r') In [{y - y'f + {x- x') 



r\2 



dS'. 



(4) 



Defining the current density variation 5J = J — J, we obtain from Eqs. and ^ that the current 
density which minimizes the functional also minimizes the functional J^' , defined as 



T'[5J] = l: I 5J{v)8Aj{v)dS+ f 5J{r)5A,{r)dS, 
2 Js Js 



Is JS 
where 6Aj is the vector potential created by 5J and 5As, = via — ^a- 



(5) 
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C. Minimization of T and magnetic energy 

In this section we demonstrate that, for the initial stage, the principle of minimization of T 
is equivalent to the magnetic energy minimization (MEM) provided that current density pene- 
trates monotonically from the surface inwards and the rate of increasing the transport current is 
proportional to the rate of increasing the applied field. 

For calculating the magnetic energy, we assume that the transport current in the strip of Fig. 
n returns through another identical one at a large distance D (D ^ a,b). In the following, we 
consider that the returning strip is centered at {x,y) = {D,0) [43|. Using the general formula for 
the magnetic energy in an infinitely long circuit W = (1/2) Jo J{r)Aj{r) + Jo J(r)Aa(r), where 
Sxy refers to the whole xy plane area, we obtain that the magnetic energy per strip, W' , is 

W' = \[ J{v)Aj{v) + [ J(r)^a(r), (6) 

ignoring constant terms that are irrelevant for MEM. W of Eq. Q is independent of the position 
of the returning strip, as long as it is placed at a enough large distance from the other strip. 

In order to compare the minimization of J^' at every time layer with the magnetic energy 



3- 



Given a physical J, we divide it 



minimization for a J in the initial stage, we do the following 
into n terms 5Ji, so that J(r) = Yli=i ^Jii^)^ being n a large number. We choose these terms as 
the actual current density increments in the initial stage corresponding to the time layer at t = ti, 
5Ji = 6J{t = ti) with ti > ti-i. If the current front penetrates monotonically from the surface 
inwards with increasing /, each 6Ji encloses a current-free core. Furthermore, if n is very high, 5Ji 
is nonzero in a thin layer only, so that the vector potential variation at time ti, 5Ai = 6Aj^i + 6Ag,^i, 
is almost uniform in the layer and 

T'[6Ji] « ^A'iili - h.i) + 1 f dJi6A^,„ (7) 



2Js 

where 6A'[ is the variation of the vector potential in the current-free core and li is the transport 
current at time ti with /j=o = 0. 

From Eqs. ©, @ and ((7)), we can separate W' as 

W' = - f J'A'jdS+ [ J'A'^dS 
2 Js Js 

+ 6Al{I-h/2) 

+ \ I 6Ji {A^ - Aa,i/2) d5, (8) 
2 Js 

where the label '1' corresponds to the first increment of current density set in the superconductor. 
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and J' = J — 5Ji, A' = A'j — 5Aj^i and A'^ = A^, — 5j4a,i. To obtain Eq. Q we used that all the 
current density created after 5Ji lies inside the current-free core of 5Ji . 
Following the same steps for all 5Ji, we find that 



UAi{2I-h-h.r) 



i=l 

+ ]- [ (5Ji(2^a-^M-^a,i-l)d5 

2 Js 



(9) 



If the final current density J is reached by increasing / and A^ monotonically from zero and the 
increase rates of / and A^ are proportional to each other, then 2A.^ — Ag^i — = 5Ai^i{2I — 

li — — /j-i). Inserting this into Eq. @ and using Eq. ((71), we obtain 

f2I-Ii-Ii.i 



w' = Y.T'm{^^j^^j^\. (10) 

i=l ^ * i 1 / 

From Eq. (jlUj) we directly deduce that when minimizing T'[6J-^ for each time layer, W' is also 
minimized, since / and /j are fixed external parameters. If the 5Ji which minimizes J-'[6Ji\ is 
unique, the J = Y^^=i ^Ji minimizing W is also unique. 

D. Calculation of current distribution 

We calculate the current distribution by minimization of J-'[&J] of Eq. ^ for each time layer 
as follows. 
As done in 



IS 



23 



4a ]. each superconducting strip is divided into N = x 2ny elements 



with dimensions a/n^ ^ b/uy] current density is assumed to be uniform in each element. In order 
to obtain a smoother current front, we allow the current density to have discrete values below Jc, 
that is, J = kJc/m with k being an integer number from 1 to a maximum value m. As discussed 



in Refs. 



48j . this reduces the discretization error in our ac loss calculations. In this paper we 
use between N = 12000 and 16000 elements and m = 20 current steps. 

Given a current /, applied vector potential Ag^ and a current density J, we calculate the current 
density variation 6 J if the current is changed into / and the applied vector potential into A^^, as 
follows. 

First, we find the element with sJ < Jc, being s = sign(/ — /), where increasing the current 
by AI = sJcoh/rixnym produces the minimum increase of J-' and we repeat the procedure until 
the current reaches I. Then, we find the elements where changing the current by AI and —A/, 
respectively, reduces the most and \ J\ does not exceed Jc; we repeat this procedure until varying 
the current in any pair of elements will increase J-' instead of lowering it. This allows the creation 
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of regions with current density opposite to /. Setting the current in this way, we ensure that the 
constrains of Eqs. ^ and Q are fulfilled. In this procedure the time do not play any role, so that 
the resulting current distribution is independent on the specific current (and field) waveform. 

The variation of J-'' , AJ^' , due to a variation of current AI in the element j can be calculated 
from Eq. © and = —^qH^x taking into account the division into elements of the tape, with 
the result 



n , 
AT'^ = Y,6IkAICjk + i:{AI)^C,j 



k=l 



- fio{H.-H,)^Ixj, (11) 

where 5/^ is the current flowing through the element k induced after the change of / and Hg^, xj is 
the X coordinate of the center of element i, and Cjk are geometrical parameters calculated in the 



appendix of 



45 



For simplicity, we consider a constant variation of / (and -f^a) between different time layers 
inside each half cycle. In this paper we use between 80 and 320 time layers. 

It can be demonstrated that the current distribution in the reverse and returning stages can 
be obtained from that one in the initial stage as long as in that stage the current front penetrates 
monotonically from the surface inwards with increasing / (and Ha) Ref. 3]. Specifically, the 
current distribution in the reverse and returning stages, Jrev and Jretj are, respectively, 

Jrev(/) = Jm(/m) " 2Jin[(/m " I)/2], (12) 
Jret(/) = -Jin(/m) + 2Jin[(/m + /)/2], (13) 



where Jin is the current distribution in the initial stage 



56| 



E. Calculation of ac loss 

The power loss per unit volume in a conductor is J • E. Then, the ac loss per unit length and 
cycle Q in the superconducting strip is 

Q = (j) dt J J{x,y;t)E{x,y;t)dxdy, (14) 

where the time integral is over one period. 

Since the electrical field inside the superconductor is in the z direction, we obtain that 

E = - A, (15) 
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where (j) is the electrical scalar potential, dz4> = d(j)/dz and A = dA/dt. The electrical field and 
the vector potential have zero components in the x and y directions, so that d(f)/dx = d(j)/dy = 0. 
Then, since for infinitely long conductors E do not depend on z, dzcf) is uniform in the whole 
cross-section. The quantity dz4> can be calculated taking one point where £" = and using Eq. 
((T3|) . obtaining that 

dz4>{t) = -A{xQ,yQ;t), (16) 

where xq and yo are the x and y coordinates at some point where E = For the critical-state 
model, E always vanishes in the positions where J = 0. Such position can be either in the current- 
free core or on the flux fronts (Sec. Ill A|l . 

Inserting Eqs. ((T31) and (fT^ into Eq. H14() and using A = IdiA and dl = Idt, we obtain 

Q = 2 dl I Jrctix,y;I) 

[diAixo,yo; I) - djA{x, y; I)]dxdy, (17) 

where the current integration is performed in the returning stage. From Eq. (|17)) . we see that Q is 
independent on the specific I{t) waveform, as long as / increases or decreases monotonically with 
time in a half cycle. From this feature it is directly deduced that the ac loss due to a sinusoidal 
current and applied field is independent on their frequency. 

The vector potential can be easily calculated from J, obtained by means of the numerical 
procedure described in Sees. Ill Bl and III Dl Then, we calculate djA at a certain time layer k from 
the numerically obtained A as 

djAix,y;I,) ^ Mx,y;Y)-Aix,y;I,^,) ^ ^^^^ 

where It is the current in the time layer k. Equation (|18|) yield much more accurate results of 
djA than using finite differences between consecutive time layers. Indeed, according to the mean 
value theorem, there must exist some current between Ik-i and Ik+i where the derivative is exactly 
the right-side part of Eq. (fTH|) . Equation (fTH|) cannot be used at the boundaries of a half cycle, 
/ = =b/in, since there diA is not continuous. Therefore, we use finite differences between k and 
k + 1 oi k and A; — 1 for / = —Im and / = /m, respectively. 

If in the initial stage the current front monotonically penetrates with increasing Jm (or i^m); 
it is possible to obtain a formula for the ac loss without derivatives in the vector potential. This 
is the case, e.g. of the pure ac transport or the magnetization by an ac external field. For this 
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situation, Jrev and J^et obey Eqs. ((T^ and (fT^ . Thanks to this feature and fohowing the same 
deduction as Carr for the pure transport situation 4^, the ac loss from Eq. (|17j) becomes 



Q = 4Jc / s{x,y) 



Ajnix,y) dxdy, 



(19) 



where and Aj^(x,y) are those corresponding to the peak values of / and -ffa and s{x,y) is a 
function giving the sign of Jrev Equation (|19)) with s = 1 corresponds to that obtained by Norris 
or the transport case [l^ and with s = x/\x\ it corresponds to the magnetic one given by Rhyner 

Since the condition of monotonic current front movement is not always valid, it is saver use Eq. 
H17() for a general treatment of simultaneous ac magnetic field and alternating transport current. 



III. RESULTS AND DISCUSSION 



In this section we present our results for the current distribution and the ac loss and we discuss 
the existing analytical approximations for low and high b/a aspect ratios. We also introduce the 
dissipation factor T = 2irQ/{iJ,ol'^), characterizing the loss behavior better than the ac loss itself. 



A. Current distribution 



In the following we present the current distribution for a rectangular strip with aspect ratio 
b/a = 0.2, although the numerical procedure gives accurate results for b/a between 0.001 and 100. 
We consider several situations of field and current. 

First, we study the case of low applied fields. As an example, we plot the current distribution 
for Im/h = 0.8 and H^/Hp = 0.08 in Fig. |2 These profiles are qualitatively similar to those for 
transport current with the difference that the current-free core is shifted to the left. In this 
situation, the current fronts monotonically penetrate from the surface inwards. Thus, the current 
distribution can be calculated using MEM (Sec. Ill C|) and the current distribution for the whole 
cycle can be constructed from that one in the initial stage using Eqs. (|12|) and ()13() . For this case, 
the ac loss can be calculated using Eq. (|19j) . so that the evaluation of E can be skipped. 

The most representative situation of combined action of ac field and ac current is that of higher 
applied fields, such as Im/Ic = 0.6 and H^/Hp = 0.72, presented in Fig. |31 The most significant 
issue is that the current distribution in the reverse stage is not always a superposition of that for the 
initial stage. Not even the returning stage is related to the reverse one. However, as can be seen in 
the figure, the current distribution for I / = 1 (and Hg^/H^ = 1) in the reverse stage corresponds 
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to that for I /Ira = —1 (and H.^/Hja = —1) for the returning one, except some numerical deviation. 
Then, the current distribution for the following reverse stage for a certain transport current / (and 
applied field Hg) is the same but with opposite sign with respect to those for the returning stage 
for transport current — / (and applied field —H^,), being current distribution periodic in time after 
the first cycle. 

The above specific case (Fig. |3| presents a current-free core, but it is not always the case for 
higher H^a or Im. For example, the current-free core is not present for Im/Ic = 0.6 and H^a/He, = 1-2 
(Fig. 0]), as well as for any case with Jm = Ic-, as shown in Fig. ISJfor /m = Ic and H^/Hp = 2. In 
addition, for all the cases with /m = Ic the electromagnetic history is erased at the end of one half 
cycle (Fig. Thus, for this current amplitude, the returning profiles correspond to the reverse 
ones with inverted sign of the current density, so that the electromagnetic behavior is simplified. 
Another issue is that the current distribution for /m = Ic have only one boundary between the 
zone of positive and negative current, while they can have two or more for lower Im (Figs. |31and 

B 

For other aspect ratios we found the same qualitative behavior as for b/a = 0.2 described above. 
As an example, in Fig. El we present the current distribution in the returning curve for b/a = 5. 
This corresponds to the situation of the same strip as for Figs. l^toElbut with applied field parallel 
to the wide direction. 



1. Comparison with analytical limits 

It is interesting to compare the sheet current density K in a thin film from Refs. [2^, where 
it is assumed that current fronts penetrate monotonically, with our results for finite thickness. We 
calculated K by integrating the current distribution over the thickness. For this situation, in Refs. 
25I . it is distinguished between the low-field high-current regime, for which all current has the 
same sign, and the high-field low-current one, when current density with both signs exist. These 
regimes appear for I/Ic > tanh(i/a/-f^c) and I / Ic < tanh(ffa/-f^c)) respectively, being He = 2Jc6/7r. 

In Fig. we present our numerical calculations of K for the initial stage together with the 
analytical results for a thin film for b/a = 0.01, Hm/Hp = 0.1 and Im/Ic = 1 (a), belonging to the 
low-field high-current regime, and Hja/Hp = 1 and Im/Ic = 1 (b), as an example for the high-field 
low-current case. As can be seen in Fig. d^a), for low applied fields all numerical results fall on 
the analytical curve within the numerical accuracy, while for higher fields, Fig. Efb), there is only 
coincidence for the profiles corresponding to low current penetration. The discrepancy for higher 
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penetration appears since the assumption of monotonic penetration of current fronts is no longer 
valid for the analytical solution. As can be seen in Fig. 2fa,b,c), for high I/Iq there is a recession 
of the zone with negative current density in favor of that one with positive current density, being 
this effect more important for higher /. The validity of the analytical solution for thin strips was 
already discussed in Refs. 2^, [2^. 

We can also compare the numerically obtained current distribution to the analytical solution 



mm- 

For strips with high 



for a slab in parallel field, for which the current fronts are planar 
b/a in the high-field low-current regime, the calculated current fronts approach to planar ones 
(Fig. El), the approximation being better for higher h/a. This behavior is in agreement to the pure 



,ppro: 



magnetic case |21l |. However, for the low-field high-current regime current fronts are similar to that 



i],Q 



and, thus, the slab 



ones for a thin strip with only transport current, which are nonplanar 
approximation is no longer valid. 

We have performed numerical simulations for very high applied fields, > and have 

shown that the current fronts approach to vertical planes for any aspect ratio, in accordance to 
the slab approximation. This is because when the applied field variation is much higher than the 
field created by the variation of J, the first term of J-' in Eq. ^ can be neglected. Since A^. is 
proportional to x, T' of the new induced current density is independent of its y location, and the 
current density profiles must be planar. 

B. Total ac loss 

First, we study the ac loss for several hja aspect ratios and their possible analytical approxi- 
mations. For this purpose, we present our results of the normalized ac loss q = 2-kQ / [iiqI'^) as a 
function of /m and constant and as a function of and constant /m in Figs. ISTallMal llUr a) 
and Figs. IHTbip^blfTnTb). respectively. Figures El IHl and EH are for aspect ratios h/a = 0.001, 100 
and 0.1, respectively. The aspect ratios of 6/a = 0.001 and 0.1 can be used to describe qualitatively 
YBCO coated conductors and Ag/Bi-2223 tapes, respectively, in perpendicular field. Aspect ratio 
b/a = 100 is representative for parallel applied field. For all figures, we consider Jm normalized to 



Ic and normalized to the full penetration field i?p, which for a rectangular strip is 
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Hp — Jc — 

TT 



2a h f a^ 

— arctan — h m 1 + 775- 
b a \ b-^ 



(20) 



The numerical error in the ac loss have been analyzed by using several numbers of elements, 
current steps and time layers, showing a variation that cannot be appreciated for the axis scale of 
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all figures below. 

From Figs. El El and EH we see that the ac loss monotonically increases with increasing either 
the current or the applied field amplitudes. For high applied field, Q increases linearly with 
for constant /m [Figs. IHtb)jntb) and liur b)]. The ac loss for the low-current limit in figures (a) and 
the low-applied- field limit in figures (b) is constant, corresponding to the pure transport and pure 
magnetic case, respectively. This qualitative behavior is consistent with experiments for YBCO 



coated conductors 



3S 



40| and Ag/Bi-2223 tapes 



31 



35 



3e 



3H. 



As expected, the loss results for zero applied field and zero transport current are the same to 
the results for pure transport and pure magnetic situations calculated using MEM in Refs. [l^ 



and 



23l |. respectively. 



1. Analytical limits for the ac loss 

In the following we study first the validity of the analytical limits for thin strips and slabs in 
parallel applied field. 

In Fig. IHlwe compare our numerical results of 2ttQ / {hqI"^) for b/a = 0.001 (line plus symbols) 
with the ac loss for an infinitely thin strip calculated by Schonborg j^^] (dash line) from the sheet 
current distribution obtained in Q and Q 

Schonborg expression for 1^ = corresponds to 
the Norris formula for a thin strip with pure transport current [l^. The thick continuous line 
plotted in Fig. |H1 separates the low-field high-current regime from the high-current low-field one. 
Sec. nil A II As can be seen in Fig. |H{a,b) the ac loss for the low-field high-current regime is well 
described by the analytical expressions for a thin strip. However, there is a significant deviation for 
the hi gh-f ield regime, increasing with increasing field or current. This is so because, as discussed 



m 



26l |. the current density formulae for thin strips are only valid for monotonic penetration 
of current fronts, which appears only for the low-field high-current regime. Fig. (21 The current 
front penetration deviates more from the monotonic case for higher field and current, so that the 
formulae for thin strips are less applicable. 

In Fig. EH we plot our numerically calculated 2Tr Q / {iiqI^) for b/a = 0.1 (line with symbols) 
together with that for a thin strip (dash line). In this figure we can see that the thin-film approx- 
imation is not valid for b/a for any case except /m = Ic and low applied field. 

It is also interesting to compare our numerical results to the formulae for the ac loss obtained 
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by Carr for a slab in parallel applied field assuming planar current fronts 24^, which in SI are 



2ttQ 
2ttQ 
27rQ 



na o 
36 



36^ 



vra , 



1+3- 



1 + 



27ra 



< /im < 1 



(1 - im)(l + «m + «m) 



(21) 
(22) 



27ra .9 (1 - 

6 '"/in, 

4vr6.2 (1 



V > 1, 



where ir 



3a {hra ^m)^ 

Ira/Ic and /im = H^J^JcO). The high-field limit of Eq. H23() is 



27rQ 



vra 



/in 



1 + 



/im > 1. 



(23) 



(24) 



In Fig. injwe plot our numerical results of 2-kQ /{fiQl^) for b/a = 100 (line with symbols) together 
with those for a slab calculated from Eqs. H2ip -()23 p (dash line). We see that the above formulae 
for slabs agree well with the numerical results for high fields and low currents, although they do 
not for low fields and high currents. In Fig. 121 we also see that for much above Hp, the Carr's 
results approach to the actual loss for any current. These features can be explained from the 
current distribution, discussed in Sec. 

The Carr formula can be also compared to numerical results for any b/a. In Figs. IHl^a) and llOf a) 
we include the high-field limit of the ac loss in a slab, Eq. (|24() . for the highest values of H^/Hp 
in those graphs (dotted lines). It can be seen that the analytical limit of Eq. (|24j) approaches to 
the numerical results for high for H^/Hp > 1 and H^/Hp > 5 for b/a = 0.001 and b/a = 0.1, 
respectively. Numerical calculations for other 6/a, like b/a = 1, also agree with Eq. H24|) for high 
applied field amplitudes. This feature can be explained as follows. For high H^a, the current fronts 
are planar, like those for a slab. Sec. IIII A II Moreover, if is high enough, the only relevant 
contribution to the electrical field, Eqs. (jlSj) and (|16|) . is from the applied vector potential for any 
aspect ratio. Then, the high-field limit for a slab must be valid for any aspect ratio. Actually, Eq. 
pi]) can be easily deduced from Eqs. (fTl|) - (|16j) assuming that A k, Aa = —fj-oxHa- 



C. Dissipation factor 



Usually ac loss under alternatin g fie 



fixing i^m or vice versa, 



27 



31 



33 



d and current have been studied as a function of Im and 
35, 36, 37, 38, 

Here we underline the significance 
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of the ac-loss dependence when increasing simultaneously /m and Hj^ with both parameters pro- 
portional to each other along the curve. This situation is found in actual ac devices, such as an 
alternating magnet. 

As explained below, for oc we can see more details of the ac loss behavior if we plot 
Q normalized to instead to I^. Actually, the quantity 2t:Q / {hqI'^ = F is proportional to the 
ac loss of a winding per the stored magnetic energy averaged during the cycle duration. Thus, F 
can be regarded as a dissipation factor. Moreover, F for only transport current is proportional to 
the imaginary part of the self-inductance, defined in j5l|, and for only applied magnetic field F is 
proportional to the imaginary part of the ac susceptibility 0, [3| • 

In Fig. ^2 we present our numerical results of F for h/a = 0.1 as a function of Im when is 
varied proportionally to /m as H^^/Hp = alja/Ic-, where a is a constant (line with symbols). This 
figure shows that for the low (and Hjo) limit, F increases proportionally to Im (or i^m)) which 
corresponds to a dependence proportional to for the ac loss. Moreover, for high a, F decreases 
with increasing with a slope in log- log scale slightly higher than —1 (and a slope around 1 for 
the ac loss), presenting a peak at a certain value of Im (or H^a). We notice that in a log- log plot 
of q against Im, the ac loss always increase with Im, appearing curves very similar to straight lines 
with a slight change in the slope. However, for F the qualitative behavior of the loss with varying 
Im and a is more evident. 

The linear dependence of the dissipation factor with Im (and Hm) for the low-field limit is 
characteristic of the critical-state model, also found for the pure transport and pure magnetic 
situations 3, l^^ • This is so since for low enough Im the current distribution is approxim ately 
parallel to the surface, as well as the magnetic field in the region with nonzero current density j2j]. 
This means that the loss factor will increase as Im (and Q as Ij^) at low levels of excitation in a 
superconductor winding of any shape or number of turns. However, for strips with very small h/a, 
such as b/a = 0.001, the linear dependence of F with Im appears only for very small Im, presenting 
for higher Im the I^ dependence typical for thin films [23 1. 

For comparison, in Fig. we also include the loss factor for only transport current (dash- 
dot line) and that for only applied field (dash lines), extracted from the tables in 

Q and Q, 

respectively, and interpolating for intermediate values of Im or Hm when needed. In a F(Im/Ic) 
representation, the curve for the pure magnetic ac loss is different for each a. We see that for low 
a, F approaches to that one for only transport for any Im and for high a, F approaches to that 
one for only magnetic field for not very high Im. 
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For a finer approximation for liigh a, we can consider tlie following loss factor 

2tt 



(25) 



where Qa^oo is an approximated ac loss as 



3/. 



(26) 



The first term of Eq. (|26|) is the ac loss with only applied magnetic field, while the second one is 
the high- field limit for an slab, Eq. (|24() . subtracting the ac loss for = 0. In Fig. ^Jwe plot 
the loss factor of E g. (|25|) for a = 2 and 5, calculated from the tables of numerically calculated 
ac susceptibility in [23| and Eq. (|25|) . We see that the approximation of Eq. (|25|) improves with 



increasing a, being almost overlapped with our numerical results for a < 5. 

For low b/a, such as 6/a = 0.001 or lower, the ac loss for a = approaches to the Norris formula 

n n 

|1RI |. if Ira is not very low |lij |. 



for thin strips [l^], if /m is not very low [13|. For the high-a limit we can obtain an 



solution of F by inserting the formula for the ac loss in a thin strip with /„ 
(dnj, obtaining 



0, Ref. 



analytical 



25|, into Eq. 



F 



2 Ha 



— I — ^ ( In cosh — - 

3 V Hr 



tanh ■ 



Hr 



(27) 



where He = 26Jc/7r. Equation ()27() is not valid for very low Hai/Hc, since the thin strip approxima- 
tion for only applied field is not valid for this range 23]. For intermediate a, F can be approximated 
from the Schonborg's formula for the ac loss in thin strips ;2^, as long as Im > tanh[7r//in/(2aJc)]) 

secnnnn 



IV. COMPARISON WITH EXPERIMENTS 



The results of our ac loss calculations presented in Figs. 



published measurements for Ag/Bi-2223 tapes 



conductors 



3£ 



31 



and [TqI qualitativelv agree with 



M M, \m M, mjm and YBCO coated 



40|. It is interesting to analyze in detail Fig. 10 of Ref. |39|. There, it is shown a 



comparison between the measured ac loss in a YBCO coated conductor and the theoretical one for 



. It can 



a thin strip in the critical state model, evaluated from the current distribution in 
be seen that the measured ac loss lay below the thin strip approximation, in agreement with our 
numerical results of Fig. |H1 As discussed in Sec. lIIIBTl this is so because the thin strip calculations 



in Refs. 



27 



39| are not valid for high applied fields. This shows that our numerical calculations 



can be used to simulate the ac loss in YBCO coated conductors. 
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In order to perform a more detailed comparison, we measured the loss factor for a commercial 
Ag/Bi-2223 tape with 37 filaments manufactured by Australian Superconductor. The sample was 
of 8cm length and 3.2 x 0.31mm cross-section with a critical current of 38A in self field at 77K. The 
superconducting core cross-section was roughly elliptical with dimensions 2a x 26 = 3.0 x 0.13mm. 
The measurements were performed at a 72Hz frequency and 77K temperature; the details of the 
experimental technique will be presented elsewhere. We present the measured results in Fig. El 
(dot line with symbols) together with numerical calculations for a rectangular strip with the same 
thickness, width and critical current (solid lines). We notice that for the theoretical curves we do 
not fit any parameter to the measured ones. In Fig. El we label the curves with the parameter 
2aH-[a/Ira instead of a in order to avoid supposing any model for performing the measurements. 
For comparison, we also include the loss factor assuming elliptical cross section for only transport 



15l | (lower dash curve) and only applied field for the highest 2aH^/I^ ISSi] (upper dash 



current 
curve) . 

From Fig. El we see that the main qualitative features of the measurements correspond to the 
behavior of a strip assuming the critical state model, except close to Ic for high 2aHi^/ 1^. This can 
be explained from the magnetic field B dependence of Jc, for which decreases with increasing 
\B\. Then, for higher //m, 1^ is lower and a normal resistive current appears in the silver for 
-^m < Ic, adding a certain contribution to F. Figure El shows that there is a better agreement 
between the measured F and that for an elliptical bar in the critical state model approximation 
than for the rectangular one, explained by the overall shape of the tape superconducting core. 
Moreover, the fact that the multifilamentary superconducting core behaves as a single solid wire 
for any suggests that the interfilamentary currents in the tape are saturated due to the high 
length of the sample 



V. CONCLUSIONS 



In this paper we have presented a systematic theoretical study of the current distribution and 
ac loss in a rectangular strip transporting an alternating transport current / in phase with an 
applied field Ha perpendicular to the current fiow. We assumed that the superconductor follows 
the critical state model with a constant Jc. With this assumption, we have developed a numerical 
procedure which takes into account the finite thickness of the strip. General features of the critical 
state model in such circumstances have been discussed. In order to do a systematic study, we have 
performed extensive numerical calculations for several aspect ratios and current and applied field 
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amplitudes, H^^ and Im respectively. Finally, we have performed measurements on Ag/Bi-2223 
tapes to be confronted with calculations. Good qualitative and quantitative agreement without 
fitting parameters have been found. 

The results for the current distribution have shown a rich phenomenology due to the metastable 
nature of the electrical currents flowing in the superconductor. For low i^a and high /, the current 
distribution is qualitatively similar to the pure transport situation. Then, J at the reverse and 
returning stages are a superposition of J in the initial one [Eqs. ()12() and However, it is not 

the same for high H^, or low / due to the nonmonotonic penetration of current fronts. In general, 
the returning stage can be deduced from that one in the first reverse stage. Only after the first 
cycle, the behavior becomes periodic. 

The ac loss Q has been accurately calculated for the thickness-to-width aspect ratios, 6/a=0.001, 
0.1 and 100, in order to qualitatively describe YBCO coated conductors and Ag/Bi-2223 tapes with 
applied fields in the transverse direction (6/a =0.001 and 0.1, respectively) and in the parallel one 
{b/a =100). Their current and applied field dependence is in accordance with published measured 
data for YBCO coated conductors Q, Q and Ag/Bi-2223 tapes Q, Q, Q, Q H, Q. We 
have shown that the ac loss behavior can be better characterized by means of the loss factor 
r = 27rQ/(^o-^m) studied as a function of /m with Hj^ proportional to /m. We have measured the 
loss factor in actual Ag/Bi-2223 tapes, obtaining a good agreement with the calculations. 

We have also presented a detailed study of the analytical limits for Q and F and their applica- 
bility. For thin samples like YBCO coated conductors, b/a < 0.001, the current profiles and the ac 
loss only approach to those for the analytical limit for thin strips [23] for the low-field high-current 
regime only. The thin film approximation is never valid for 6/a ~ 10, such as for Ag/Bi-2223 tapes. 
We have also studied the slab limit, obtaining that for the situation of parallel field, b/a = 100, 
the slab approximation is not valid for the transport-like regime (low-ffm arid high-/m)- However, 
the high-field limit for the slab approximation can be used for any aspect ratio provided that 
is high enough. 
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FIG. 1: Sketch of the tape cross-section and division into elements for the calculations. The transport 
current / and the applied field Ha are directed in the positive z and y directions, respectively. 

FIG. 2: Current distribution in the initial stage for the low-field and high-current regime. Specific parameters 
are hja = 0.2, H^/Hp = 0.08, I^/Ic = 0.8 and 7//^ = H^/Hp = 0.2 (a), 0.6 (b) and 1 (c). Local current 
density is + Jc for the black region and zero for the white one. 

FIG. 3: Current distribution for b/a = 0.2, H^/Hp = 0.72, and Im/Ic = 0.6 at several instants of the ac 
cycle. Figures (a,b,c) are for the initial stage with I/Im = 0.2, 0.6, and 1, respectively, figures (d,e,f) are for 
the reverse stage with I/Im = 0.6, —0.2 and —1, respectively, and figures (g,h,i) are for the returning stage 
with I/Im = —1, 0.2, and 1, respectively. In the black regions the local curent density is +Jc, in the light 
gray zones it is — Jc, and in the white ones it is zero. 
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FIG. 4: Current distribution for b/a ~ 0.2, Hn^/Hp = 1.2 and Im/Ic — 0.6 at several instants of the ac 
cycle. Figures (a,b,c) are for the initial stage with 1/1^ = 0.2, 0.6, and 1, respectively, and figures (d,e,f) 
are for the reverse stage with ///,„ — 0.6, —0.2 and —1, respectively. In the black regions the local curent 
density is + Jc and it is —Jc in the gray zones. 

FIG. 5: Current distribution at the reverse stage for b/a — 0.2, H^/Hp ~ 2, Im/Ic ~ 1 and I/Im ~ 0.6 (a), 
—0.2 (b), and —1 (c), respectively. In the black regions the local curent density is +Jc and it is — Jc in the 
gray zones. 



FIG. 6: Current distribution at the reverse stage for b/a — 5, Hm/Hp = 0.72, Im/Ic = 0.6 and I/Im ~ 0.6 
(a), —0.2 (b), and —1 (c), respectively. In the black regions the local curent density is +Jc and it is — Jc in 
the gray zones. 

FIG. 7: Sheed current density K in the initial stage as a function of x for b/a = 0.01, Im ~ Ic, and 
Hm/ IIn = 0.1 (a) and 1 (b) at several instantaneous / (and iJa)- Lines are for the thin strip limit from 
Refs. [2^ |2a| and symbols are for our numerical calculations. For the numerical results, K is the integral 
of J over the sample thickess. 



FIG. 8: Normalized ac loss 2ttQ / {i.iqI^) for b/a = 0.001 as a fuction of Im/Ic for several Hm/Hp (a) and 
as a fuction of Hm/Hp for several Im/Ic (b). Lines with symbols are for our numerically calulated results, 
dash lines are for the thin strip limit from 12 7|. the thick solid line separates the low-field and high-current 
regime from the high-field one in a thin strip l27|, and dotted lines (for Hm/Hp=l and 2) correspond to the 
high-field limit for slabs [Eq. if^ ]. 



FIG. 9: Normalized ac loss 27rQ/(^o4^) for b/a = 100 as a function of Im/Ic (a) and as a fuction of 
Hm/Hp (b). Solid lines with symbols are our numerically calulated results and dash lines are for the slab 
approximation from Eqs. IHJ-lISl) Q- 



FIG. 10: Normalized ac loss 2ttQ/{ij.oI^) for b/a = 0.1 as a fuction of Im/Ic (a) and as a fuction of Hm/Hp 
(h) . Lines with symbols are for our numerically calulated results, dash lines are for the thin strip limit from 
[27|, and dotted lines (for iJ,„/iJp=2 and 5) correspond to the high-field limit for slabs [Eq. H24ll ]. 



FIG. 11: Loss factor F = 2ttQ/ (/ioJ^) as a function of Im/Ic with Hm proportional to Im as Hm/Hp ~ alm/Ic 
for several a. Solid lines with symbols are for our numerical calculations, the dash-dot line corresponds to 
only transport current, the dash ones are for only applied magnetic field, and dotted lines correspond to the 
high-a approximation from Eqs. H25(l and l|26(l . 
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FIG. 12: Calculated loss factor T together with experimental data from a commercial Ag/Bi-2223 tape. 
Dot lines with symbols are for measurements, solid lines correspond to numerical calculations assuming a 
rectangular cross-section, and the two dash ones are for elliptical cross-section for i/m = (lower curve) and 
— (top curve) from Refs. Q and 0, 

respectively. 
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